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SUMMARY 

I.  Mode  Locking  and  P>'!se  Selection 

In  the  previous  semiannual  report,  work  of  Y.  Manichaikul  was 
described  in  which  he  suet  essfully  switched  out  a  single  optical  pulse 
from  a  mode -locked  C02  transversely  excited  atmospheric  pressure  (TEA) 
laser.  The  peak  pulse  power  vms  of  the  order  of  10  kW. 

A  rew  system  has  been  designed  to  produce  short  L.ser  pulses  in 
order  to  improve  the  achievable  peak  powers  and  also  to  enable  us  to 
obtain  two  optical  pulses  of  roughly  equal  intensity  separated  by  a  time 
of  the  order  of  20  ns.  Improvements  in  the  peak  powers  were  necessary 
to  enable  saturation  of  a  TEA  amplifier  with  the  pulse  in  order  to  study 
nonlinear  amplification  in  such  a  system;  the  selection  of  two  (or  more) 
successive  pulses  was  desirable,  in  order  to  ascertain  the  relaxation 
times  responsible  for  the  recovery  of  the  upper  laser  level  through  inter¬ 
change  of  energy  with  other  vibrational  levels  of  file  same  vibrational 
(asymmetric  stretching)  mode.  The  system  has  been  operated  succes 
fully,  giving  pulses  of  less  than  2  ns  duration  and  501  kW  peak  intensity. 

The  system  is  described  in  detail  in  Quarterly  Progress  Report  No.  110, 


Research  Laboratory  of  Electronics,  July  15,  1973;  a  revised  version 
meantime  appears  as  Appendix  I  of  this  report.  We  have  now  made 
extensive  amplifier  measurements  which  we  are  analyzing  and  will 
report  subsequently. 

In  addition  to  forced  mode  locking  by  an  acoustically  driven  crystal, 
mode  locking  by  a  saturable  absorber  (a  hot  CC>2  cell)  promises  to  be 
a  convenient  means  to  produce  short  mode-locked  dulses.  Such  work 
on  conventional  TEA  lasers  has  been  reported  in  the  literature.  The 
scheme  is  particularly  attractive  for  the  mode  locking  of  lasers  oper¬ 
ating  at  10  atm.  Such  a  laser  has  recently  been  operated  successfully 
in  this  research  group  by  J.  L.  Miller.  A  gas  cell  has  the  obvious 
advantage  over  a  crystal  in  that  it  is  less  susceptible  to  optical  damage, 
which  is  a  more  critical  requirement  at  the  high  intensities,  achieved 
with  high-pressure  lasers.  The  theory  of  mode  locking  by  saturable 
absorbers,  although  widely  discussed  in  the  literature,  ir  less  well 
understood  than  the  theory  of  forced  mode  locking.  In  an  attempt  to 
predict  the  operating  requirements  of  the  hot  C02  cell,  a  theoretical 
analysis  has  been  initiated  of  mode  locking  by  saturable  absorbers. 

Thus  far,  we  have  analysed  the  case  of  a  weakly  saturated  implifying 
medium  and  a  weakly  saturated  absorber,  the  latter  with  response 
times  fast  compared  with  mode -locked  pulse  time  and  amplifier  re¬ 
sponse  time.  The  theory  is  presented  in  Appendix  II,  a  preprint  from 
Quarterly  Progress  Report  No.  110,  Research  Laboratory  of  Elec¬ 
tronics,  July  15,  1973. 
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II.  Closed-Form  Analysis  of  Electron  Distribution  and  Pumping  Rates 


The  analysis  of  the  electron  distribution  in  a  discharge  excited  mole¬ 
cular  laser  mentioned  in  our  previous  semiannual  report  has  been  written 
up  in  a  paper  submitted  to  the  Journal  of  Applied  Physics;  a  copy  of  this 
paper  constitutes  Appendix  III.  Stated  briefly,  the  objective  of  the  work 
is  to  find  simple  analytic  expressions  for  the  electron  distribution  in 
the  laser  discharge,  including  the  effect  of  lasing.  Examples  treated 
in  the  paper  include,  for  an  E-beam  C02  laser,  the  evaluation  of  the  V-I 
characteristic  of  the  sustainer  field,  the  gain,  the  efficiency  and  the 
pumping  rates  of  the  molecular  levels  with  and  without  lasing. 


APPENDIX  I 


Revised  Version  of  Preprint  from  Quarterly  Progress  Report 

No.  110,  Research  Laboratory  of  Electronics,  M.  I.  T. ,  July  15, 

1973,  pp.  118-121 

GENERATION  AND  AMPLIFICATION  OF  HIGH-INTENSITY 
NANOSECOND  TEA  C02  LASER  PULSES 

Y.  Manichaikul 

We  have  reported  previously  on  generation  of  short  nanosecond 
pulses  from  a  pin  resistor  TEA  C02  laser  by  way  of  mode  locking  and 
cavity  dumping.1  Single  pulses  4  ns  wide  (full  width  at  half  maximum, 
FWHM)  with  peak  power  ~10  kW  were  produced.  Problems  were  en¬ 
countered  with  the  system  and  a  new  one  was  built.  The  problems,  and 
the  changes  that  have  been  made  in  the  new  system,  are  as  follows. 

(i)  The  peak  power  obtained  from  the  previous  system  was  too  low 
for  some  experiments.  For  example,  to  saturate  a  TEA  C02  laser 
amplifier  requires  a  peak  intensity  greater  than  100  kW/cm2,  In  our 
new  system  we  use  a  3-electrode  TEA  C02  discharge  tube  which 

has  a  higher  gain  than  the  pin-resistor  TEA  CO?  discharge  tube.2  With 
this  improved  system  we  have  obtained  mode-locked  pulses  of  peak  powers 
in  excess  of  500  kW. 

(ii)  A  germanium  acousto-optic  modulator  at  Brewster  angle  is 
now  being  used.  This  eliminate?,  the  power-density  limitation  that 
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arose  in  the  old  system,  which  used  an  antireflection-coated  germanium 
modulator. 

(iii)  An  antireflection-coated  GaAs  crystal  is  now  placed  outside, 
rather  than  inside,  the  optical  cavity,  which  makes  alignment  of  the 
cavity  less  critical. 

(iv)  In  the  old  system  the  high-voltage  supply  to  the  cavity-dumping 

GaAs  crystal  was  falsely  triggered  by  electrical  noise  generated  by  the 

3 

laser  discharge  tube.  Ir  the  new  system  a  laser-induced  spark  gap  is 
used  to  eliminate  this  problem. 

(v)  Previously,  continuous  RF  power  was  supplied  to  the  germanium 
mode-locking  crystal.  This  caused  the  germanium  crystal  to  heat  up, 
thereby  changing  its  acoustic  resonance  frequencies.  We  are  now  using 
pulsed  RF  power  to  minimize  this  problem. 

Experimental  Arrangement 

Figure  1  is  a  schematic  diagram  of  the  experimental  arrangement. 
With  a  dc  voltage  cf  5  kV  across  the  GaAs  electro-optic  modulator 
a  fraction  of  the  energy  from  the  beam  can  be  switched  out  at  the 
Brewster-angle  germanium  plate,  Ge  P.  This  fraction  of  the  beam  was 
guided  by  two  flat  mirrors  to  a  gold-coated  mirror  with  2-m  radius  of 
curvature,  which  focuses  the  beam  inside  the  3 -electrode  laser  ampli¬ 
fier.  The  tube  is  operated  between  10  Torr  and  400  Torr.  A  frac¬ 
tion  of  the  beam,  before  going  into  the  amplifier,  is  reflected  out 
from  a  be.  in  splitter  and  detected  by  a  copper-doped  liquid-helium-cooled 
detector  with  a  rise  time  of  <1  ns.  The  fraction  of  the  beam  that  has 
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passed  through  the  amplifier  is  detected  by  a  gold-doped  liquid-nitrogen- 
ccoled  detector  with  a  rise  time  of  ~1  rs. 

Experiments 

a.  Generation  of  Nanosecond  TEA  CO^  Laser  Pulses 

The  laser  cavity  has  an  optical  length  of  approximately  1.88  m.  To 
achieve  forced  mode  locking,  RF  driving  .  ower  over  time  interval'-  of 
2.  5  ms  duration  with  4  W  peak  power  was  supplied  to  the  germanium 
acousto-optic  modulator.  The  3 -electrode  discharge  tube  inside  the  cav¬ 
ity  is  set  to  trigger  af  2  ms  after  the  RF  power  is  on.  We  have  obtained 
mode-locked  pulses  <2  ns  v/ide  (FWHM)  with  a  peak  power  of  ~500  kW. 

The  switching  out  of  individual  pulses  is  accomplished  by  using  the 
laser-induced  spark  gap,  which  is  normally  filled  with  prepurified  nitro¬ 
gen  at  100  psi.  The  coaxial  cable  was  charged  up  to  15  kV.  We  can 
vary  the  temporal  triggering  of  this  gap  by  altering  the  nitrogen  pres¬ 
sure,  the  distance  between  the  two  electrodes  of  the  spark  gap,  or  both. 
When  the  spark  gap  is  triggered,  a  square  voltage  pulse  of  7.  5  kV,  twice 
the  length  of  cable,  is  produced.  This  pulse  is  supplied  to  the  GaAs  crys¬ 
tal,  which  rotat<  o  the  polarization  of  the  desired  number  of  mode-loci  ed 
pulses.  These  mode  -locked  pulses  are  then  reflected  out  of  the  train  at 
the  Brewster-angle  germanium  plate. 

b.  Amplification  of  Nanosecond  TEA  CO^  Laser  Pulses 

Three  pulses  from  a  train  of  m'>de-locked  pulses  were  switched 
out  and  guided  through  an  amplifier  tube.  These  pulses  were  measured 
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before  and  after  they  propagated  through  the  laser  amplifier.  Figure  2 
shows  input  and  output  signals  detected  by  copper-doped  germanium  detec¬ 
tors.  The  output  signal  was  artificially  delayed  100  ns  by  use  o'"  a  coaxial 
delay  line.  By  using  the  add -mode  on  a  Tektronix  7904  oscilloscope,  we 
were  able  to  display  both  signals  on  a  single  trace. 

In  cur  preliminary  studies  we  have  found  that  the  first  pulse  of  our 
three  laser  pulses  is  of  sufficient  power  density  to  saturate  the  laser 
amplifier.  Further  work  on  this  aspect  of  the  experiment  will  be  described 
in  a  future  report. 
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APPENDIX  II 


Preprint  from  Quarterly  Progress  Report  No.  110,  Research 
Laboratory  of  Electronics,  M.  I.  T.  ,  July  1973,  pp.  21-23 


MODEL  OF  MODE  LOCKING  WITH  SATURABLE  ABSORBER 
C.  P.  Ausschnitt,  H.  A.  Ilaus 

Of  the  many  studies  of  mode  locking  using  saturable  absorbers,  the 
work  of  Kuizenga  and  Siegman1  on  forced  mode  locking  has  the  greatest 
analytic  simplicity.  In  this  report  our  objective  is  to  obtain  simple 
expressions  for  the  intensity  and  the  pulse  width  of  passively  mode -locked 
pulses  and  to  choose  a  model  simple  enough  to  make  closed -form  analysis 
possible.  We  assume,  with  Kuizenga  and  Siegman,1  that  the  gain  of  the 
laser  medium  changes  negligibly  through  passage  of  a  single  pulse.  We 
«  further  assume  that  the  response  time  of  the  saturable  absorber  is  fast 

enough  so  that  the  population  difference  between  the  lower  and  upper  levels 
in  the  saturable  absorber  can  follow  the  pulse  envelope  instantaneously. 
Finally,  we  assume  that  the  population  perturbation  in  the  saturable 
absorber  is  small  enough  that  its  rate  equation  may  be  linearized. 

Our  analysis  is  analogous  to  Kuizenga  and  Siegman' s  and,  in  particu¬ 
lar,  we  require  that  a  Gaussian  pulse  reproduce  in  magnitude  and  width 
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after  passage  through  the  amplifier  and  the  absorber  and  after  reflection 
from  the  output  mirror. 

Adopting  the  notation  of  Ku^enga  and  Siegman  and  following  their 
analysis  of  the  passage  of  a  Gaussian  pulse  through  the  amplifying 
medium,  we  obtain  the  electric  field  at  f,he  amplifier  exit. 


E  = 


EG 

o 

4  si y  A 


exp[-(t-B)2/4A]. 


(1) 


When  this  puJ^e  passes  through  the  saturable  absorber,  it  affects  the 
population  difference  between  the  lower  and  upper  levels  of  the  absorber 
n  according  to  the  rate  equation 


ii  + 


-n 
"V  i 


*; 


(2) 


All  symbols  in  this  equation  refer  to  the  parameters  of  the  r.aturable 
absorber.  Here  n  is  the  steady-state  value  of  n. 

Because  the  main  pulse-shaping  occurs  near  the  maximum  of  an 
optical  pulse,  we  expand  the  E -field  of  Eq.  1  around  the  instant  at  which 
E  reaches  a  maximum  and  the  absorption  reaches  a  minimum.  In  this 
way,  we  obtain  for  the  population  difference  in  the  absorber 


n  e  n 


e 


e2g2 

^-[l-2(t-B)2/4A]  •. 


(3) 


The  signal  passing  through  the  absorber  is  multiplied  by 

exp(-a£)  =  exp(-Sn),  (4) 
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where  S  is  a  parameter  characteristic  of  the  relaxation  Mme  anc  matrix 

elements  of  the  absorber.  At  this  point  cur  analysis  connects  with  the 

Kuizenga -Siegman  analysis,  since  the  modeling  of  the  absorber  makes 

it  perform  a  functio-  identical  to  the  function  performed  by  the  modulato 
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used  in  forced  mode  locking.  The  parameter  Kuizenga  anc. 

Siegman1  is  replaced  by 


T  ?  4  ' 

Is  l6yA^ 


where  I  is  now  the  saturation  intensity  in  the  absorber.  Using  the 
s 

results  of  Kuizenga  and  Siegman,  from  the  self-consistency  requirement 
on  the  pulse  width  we  find 


1  x  1  G 
Y  =  4A  +  I  2  ' 

4A  s  l6yA^ 


The  self-consistency  requirement  on  the  pulse  amplitude  gives 


exp(-Sne)  =  1, 

2sTyA 


The  parameter  R  is  the  reflectance  of  the  output  mirror.  Solving  these 
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two  equations  for  y  and  the  inter  sity  I,  and  using  A  =  l/(4y)+4g/Ao>  , 
we  obtain 


**  . t  2,-,  2  ,  i  Ad 

y  =  (L  C  -!)T^ 


~  4  In  (R/L)(L2G2-1)L2, 


t  V  .— *• 


where  we  have  defined  a  loss  parameter 


L  h  R  exp(-Sne).  (io) 

j  he  pulse-width  parameter  of  Eq.  10  can  be  rewritten  as  a  pulse  duration 
time  t  ,  which  becomes 


Tp  =  S\l  2  In  2 


/Z 


an 


The  parameters  in  Eq.  1 1  are  identical  to  those  used  by  Kuizenga  and 
Siegman,  with  the  exception  of  the  loss  parameter  L,  which  now  involves 
both  the  attenuation  of  the  saturable  absorber  and  the  reflectance  of  the 
mirror,  and  Ig  which  is  the  saturation  intensity  of  the  absorber. 

These  results  may  be  compared  with  published  experiments  on  a 
TEA  CC>2  laser  that  is  mode -locked  via  a  saturaMe  absorber  consisting 
of  a  heated  C02  cell  at  a  pressure  higher  than  tht  amplifier  pressure.*” 
The  assumptions  here  are  generally  consistent  with  these  experimental 
conditions.  The  observed  pulse  length  was  4  ns.  If  we  introduce  parame¬ 
ters  that  seem  to  match  experimental  results,  we  predict  a  pulse  length 
~1.  5  ns.  We  believe  that  better  agreement  with  experiment  can  be 
obtained  if  we  generalize  the  analysis  presented  hare  to  transient  mode 
locking,  by  using  an  analysis  analogous  to  Feldmun's  of  transient  forced 
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mode  locking. 
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ELECTRON  DISTRIBUTIONS  IN  GAS  LASERS 
W.  P.  Allis*  and  K.  A.  Haus** 

Abstract 

Simplified  models  of  collision  processes  are  introduced 
which  make  possible  closed  form  analyses  of  electron  distri¬ 
butions  in  molecular  g*ses.  In  particular,  models  are  con¬ 
sidered  that  approximate  the  collision  processes  in  a  diatomic 
gas  and  in  a  gas  mixture  resembling  the  He:N2:C02  laser. 
Expressions  arc  obtained  for  the  rates  of  excitation  of  the 
vibrational  modes  which  include  the  effect  of  superelastic 
collisions.  The  results  are  illustrated  by  evaluating  the 
V-I  characteristic  and  the  efficiency  of  an  E-bean  C02  las?r. 

Introduction 

In  an  effort  to  predict  the  gain  and  lasing  efficiency 
of  '-.w  and  CC2  lasers,  computer  studies  have  bee  1  carried  out 
to  obtain  the  electron  distribution  and  from  it  the  pumping 
raies  in  typical  gas  laser  mixtures  [1-5] .  With  such  studies 
it  is  possible  to  predict  observed  gains  and  efficiencies  as 
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functions  of  excitation  and  gas  mixtures. 

Even  though  such  computer  studies  are  now  available 
numerical  predictions  can  be  made  for  any  specific  choice  of 
parameters,  it  is  still  of  interest  to  find  models  of 
lasers  which  permit  closed  form  solutions  and  maka  possible 
easier  identification  of  the  influence  of  the  various  para¬ 
meters.  Here  we  present  a  model  of  interaction  of  electrons 
with  diatomic  and  polyatomic  molecules  which  permits  closed 
form  evaluation  of  the  electron  distribution.  The  population 
distributions  of  the  molecul.  s  over  their  vibrational  levels 
affect  the  electron  distribution  and  vice  versa.  These  ef¬ 
fects  are  brought  out  clearly  in  the  set  of  equations  ex¬ 
pressing  the  excitation  rates  in  terms  of  electric  field  and 
the  molecular  population  distribution.  The  .results  obtained 
for  a  diatomic  laser  are  compared  with  computer  results  of 
Nighan  [1]  and  we  find  reasonable  agreement.  A  two -temperature 
model  is  developed  for  a  mixture  of  C02  and  nitrogen  and  is 
applied  to  obtain  expressions  for  the  V-I  characteristic  and 
efficiency  for  an  E-beam  C02  laser . 

2.  The_Electron_Ga:j 

Electrons  driven  by  a  field  E  through  a  gas  background 
tend  to  maintain  a  spherically  symmetric  velocity  distribution 
as  long  as  the  electric  field  is  not  excessive.  This  is  the 
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case  for  the  values  of  E/N  (where  N  is  the  number  density 
of  the  gas)  employed  in  gas  laser  discharges.  We  describe 
the  electron  distribution  by  the  distribution  function  f (v) 
where  v  is  the  vector  velocity .  In  the  usual  analysis  of 
f ( v)  one  separates  f (v)  into  a  spherically  symmetric  part 
f0(v)  (where  v  is  the  racial  coordinate  In  velocity  .ipace, 
v-  |vj)  and  a  part  .^(E-v'  which  is  the  first  higher  order 
expansion  term  of  f(v)  In  a  scrie?  of  spherical  harmonics. 

The  nteady  state  Boltmann  equation  leads  to  two  coupled  first 
order  differential  equations  in  v  for  fg  and  f^.  Elimination 
of  f1  and  disregard  of  the  elastic  collision  losses  leads  to 
the  following  second  order  differential  equation  for  fg(v)[l,  6]: 


4ir 


3v 


v 

r 


e  E2  df0 

m  v“  Sv” 
*  c 


( 


W  inel 


where  v  is  the  frequency  of  momentum  transfer  collisions.  Tt.2 

c 

term  on  the  left  hand  side  represents  the  rate  of  entry  of 
elec  hr  on  s  into  the  velocity  ra.ica  v,  v  +  dv  produced  by  the 
driving  electric  field.  The  rate  of  entry  of  electrons  is 
equal  to  the  rate  of  exit  caused  by  the  inelastic  collisions, 
the  term  on  the  right  hand  side  of  the  equation.  The  distri¬ 
bution  function  fQ(v)  is  normalized 

/4tt  v2  fg  (v)  dv  »  1  W 

It  it  convenient  to  use  as  independent  variable  che  energy  in 
electron  volts: 
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(3) 


u  *  IT- 

Further,  it  is  convenient  to  renormalize  fg(u)  so  that 

/fQ(u)  /u  du  *  1  , 

a  normalization  which  avoids  factors  of  2*  and  2e/m  in  the 
equations.  We  shall  henceforth  omit  the  subscript  0.  Ncce 
that 

2e  V2 

f(u)  -  2*  (--)  f(v)  . 
m 

Wc  do  not  introduce  a  new  symbol  for  f(u),  with  the  under¬ 
standing  v.hat  the  variable  used  will  identify  which  f  is 
being  employed.  Renormalized,  ana  with  the  change  of  vari¬ 
able,  Eq.  (1)  becomes  * 


I_{2e(E_2v  u3/2  df} 

3u  'T  m  c  3uJ 


-  s 


inel 


(4) 


where  the  symbol  Sinel  is  to  imply  a  source  of  electrons  per 
unit  energy  range  as  caused  by  inelastic  collisions.  The  term 
in  the  wavy  brackets  is  conveniently  identified  with  a  function 


_  2  e  ,E  .2 

GE  "  T  n  <~)  • 


«,/2§f 


(5) 


This  function  has  the  interpretation  of  "current"  in  (one 
dimensional)  energy  space.  Indeed,  it  obe} s  the  one  dimen¬ 
sional  continuity  equation: 
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(6) 


dGE 

3u  "  sinel 

where  is  a  source  rate  per  unit  energy.  The  analogy 

with  a  "current"  can  be  carried  quite  far.  The  divergence 

of  a  current  density  is  equal  to  a  source  of  particles.  In 

the  present  case,  the  role  of  source  is  placed  by  S. 

1  2  tr»el 

which  represents  the  introduction  of  electrons  per  unit 
energy  range  via  the  inelastic  collisions.  This  situation 
illustrated  schematically  in  Fig.  1.  One  may  compare  G„ 
tu  a  conduction  current  density  J  driven  bv  an  electric  field 
through  a  medium  of  conductivity  a  .  The  "drive"  in  the  pre- 
sent  case  is  not  E,  but  E  .  The  quantity  analogous  to  a 
conductivity  is 

_2  e  u3/2  df 

7  m  v*  3u 
c 

This  "conductivity"  is,  in  general,  positive.  However,  when 
>  0,  then  instabilities  may  occur  in  energy  space,  very 
much  like  the  instability  in  real  space  associated  with  a  nega¬ 
tive  conductance. 

Next  consider  inelastic  collisions  leading  to  electro-  ic 
excitation.  Suppose  that  electrons  excite  electronic  levels 
of  the  atoms  and/or  molecules  when  they  possess  a  sharply  de¬ 
fined  energy  ux<  Electrons  are  being  taken  out  of  the  energy 
range  ux,  ux  +  du  and  we  assume  that  they  are  returned  with  a 
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smaller  energy  at  ux  ■  ux  "  Vx'  Vx  *s  the  eneray  of  t*ie 
excited  transition.  A  source  function  may  be  defined  for 
this  process ,  S  ,  which  gives  the  rate  of  entry  of  elec- 
trons  due  to  the  electronic  excitation  process  •  If  we  fur¬ 
ther  assume  that  all  electronic  excitations  take  place  at 
ux,  this  becomes  a  barrier  beyond  which  electrons  cannot  go , 
and  f  ■  0  for  u  t  ux.  ~ 

Sx(u)  -  -!vx  5  <u  -  u^  “  vx  «  (u  -  ux  +  V8)  ]  (7) 

Here,  v  is  the  electronic  excitation  frequency  and  6  is 
X 

the  Dirac  delta  function. 

In  CO,  the  excitation  of  vibrational  levels  by  the  elec¬ 
trons  occurs  through  the  formation  of  a  negative  ion  state. 
The  excitation  of  the  state  requires  1.8  electron  volts  of 
energy.  Upon  excitation,  the  electrons  lose  an  energy tro»a 
where  hoj  is  the  vibrational  energy  interval,  if  the  excita- 

Cl 

tion  carries  the  molecule  ^rriri  the  0  state  to  the  1  state. 

The  energy  loss  is  approximately  2hu?  if  the  molecule  is  ex¬ 
cited  from  the  0  to  the  2  state,  etc.  We  shall  model  the 
vibrational  excitation  process  fcy  assuming  that  all  elec¬ 
trons  lose  an  energy  Va  “  ITwa  upon  an  exciting  collision  with 
a  molecule,  when  they  reach  an  energy  ua  +  V&  (where  ua  =  1-8 
for  CO) ;  V  will  be  assumed  small  enough  when  evaluating  the 
differential  equation  for  the  distribution  function,  so  that 
the  source  function  Sa  associated  with  molecular  excitation, 


Sa(u)  ■  -tva  Mu  -  ua)  -  va  6  (u  -  ua  +  va)l 

(8) 

may  be  approximated  by  a  "doublet" ,  a  derivative  of  a  delta 
function,  (v„  •,  Va  0)  .  Note  that  v.  represents  the  net 

vibrational  "excitation"  frequency  due  to  both  excitation 
and  de-excitation.  Now,  returning  to  Eq.  (4)  we  introduce 
explicitly  the  two  source  functions  Sa  and  Sx. 


dGE  _  d 
cTu*“  “  3u 


,2  e  ,E  .  2  3/2  df, 

<3  m  ‘sH  ve  u  35> 


-  Sx(u)  +  Sa(u) 


(9) 


A  first  integral  c  ;  the  equation  gives 


C’E  -  4  I  <5->  vc  uV2  H"  l  (sx‘  »  +  Sa<u>)du 
c  o 

(10) 

The  integral  of  the  source  function  is  conveniently  identified 
with  a  symbol  of  its  own 

£  Sx(u)  du  =  -Gx(u)  dD 


/  Sa(u)  du  =  -Ga(u)  (12) 

o 

These  G  functions  themselves  may  be  identified  as  currents  in 
energy  space  whose  divergence  is  responsible  for  the  source 
of  electrons.  Indeed,  the  inelastic  collisions  also  shuttle 
the  electrons  through  energy  space.  A  net  arrival  rate  of 
electrons  per  unit  energy  is  then  caused  by  the  negative 
divergence  of  the  "currents"  G&  and  Gx  which  are  produced  by 
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the  inelastic  collisions.  Indeed,  Eqs  .  (11)  and  (12) 
may  also  be  read  as 


dGa  .  -S 

an-  s. 

One  may  now  interpret  Eq.  (10)  as  requiring  that,  in  the 
steady  state,  all  currents  in  energy  space  must  add  up  to 
zero. 

GP  +  G,  +  G  *  0  (13) 

u  A  X 

This  is  shown  schematically  in  Fig.  2.  Sinel  may  be  thought 

to  be  caused  by  (G&  +  Gx)  which  then  in  turn  causes 

It  is  convenient  to  think  in  terms  of  the  functions  G& 

and  Gx  rather  than  in  terms  of  and  Sx.  The  former  are 

shown  in  Fig.  3.  G_  is  approximated  by  a  unit  impulse  (delta) 

function  whose  area  (content)  may  be  obtained  through  simple 

physical  reasoning.  The  energy  supplied  by  the  electrons 

through  the  collisions  represented  by  S&  is 

u  u  dG  u  u 

/  u  S  du  ■  /  u  du  »  u  G  |  -  /  Gadu 
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am-m 


. . —  -  " 


/G.du  is  the  content  of  the  delta  function  G  . 

A  A 


Prom  the  above  we  may  interpret 


-  /  G-du  as  the  energy  vava  lost  by  the  electrons 


■ 


Equation  (13)  may  be  rewritten  as  a  differential  equation 


for  f(u) 


-4  -  (-  )2  v  u3/2  4“  “  "(G  +  G  ) 

3  m  c  du  '  x  a1 


(15) 


where  G  +  G  is  a  known  function  consisting  of  the  rectangle 
of  height  vx  and  the  delta  function  of  area  vaVft  shown  in 
Figure  3  .  In  order  to  integrate  it  is  necessary  to  know 
the  collision  frequency  vc  as  a  function  of  u.  For  simpli¬ 
city  we  shall  assume  it  *to  be  constant  and  then  the  integral 
is  easy  and  yields: 


f  (u)  "  vl  2S7 


—  »  ...  { — —  -  — i — }  ;  u  <  u  <  u 
c  *5  V5  /u L  a  x 


(16) 


3 _  rl  1  .  a  a _ <  «  /  u  <  u. 

”c  ^  vx  2u,3/2  ’  *  "  * 


where 


Uj  -  i  £  (S_,2 

d  7  m  '  v 

c 


(17) 


is  the  drift  velocity  of  the  electrons  expressed  in  electron 

2 

volts.  Note  that  u^  is  proportional  to  the  ratio  (E/n) 

The  normalization  condition  (1)  when  utilized  gives 


Vx  +  Va 


v  2u, 


(18) 
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This  is  the  law  of  energy  conservation  per  particle.  On  the 

left  hand  side  appears  the  rate  of  anergy  transfer  to  the 

electronic  and  vibrational  excitations,  on  the  right  hand 

2 

side  is  the  power  supplied  by  the  electric  field  (u^  «  E  )  . 
One  may  combine  (16)  and  (18)  to  eliminate  *a  and  obtain 


f  (u)  -  ^  ux  ’  u  >  ua 

_  ,V  x  .  3  (1  _  1  . 

(5Vo’  \  ^ 


2u 


'  (1  -  u  >  o  >  u  -  V  (19) 

2u  ,v  a  xx 

i  a  c 


Note  that  vxvx/2U(jV^  is  the  fraction  of  the  power  absorbed 
which  goes  into  electronic  excitation  and  must  be  less  than 
unity.  The  larger  this  parameter,  the  larger  f(u)  in  the 
range  um  <  u  <  u  .  Figure  4  shows  plots  of  f  (u)  for  dif- 

a  X 

.  v  V 

ferent  values  of  the  excitation  parameter  Exc  =  -x  X--  . 

2u.v 

The  following  values  have  been  chosen:  ug  *  2  eV,  u  =  io  ev 

v  ■  .2  eV,  V  *=10.  The  entries  E/N  will  be  explained  later, 
a  x 

One  consequence  of  the  simplifying  assumptions  is  immediately 


evident.  The  electron  distribution  is  pim.ed  down  at  u  =  ux 
and  is  zero  beyond.  (If  we  had  allowed  for  excitation  levels 
above  u  ,  this  would  not  have  happened)  .  Therefore  our 
model  is  limited  to  F/N  values  for  which  most  of  the  electrons 
have  energies  less  than  ux  and  does  not  give  information  on 
electron  densities  at  u  >  uy.  To  obtain  snch  information, 
for  example,  to  compute  ionization  rates,  one  must  do  addi¬ 
tional  modeling  as  shown  in  section  .  Another  interesting 
consequence  of  the  model  is  that  the  electron  energy  parameter 
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as  defined  by  Nighan 


u r  =  j  /  u3'2  f (u)  du 


lies  between  the  two  extremes: 


lim  u  '■* 
E+0 

and  lim 


In  particular,  in  the  limit  of  zero  field  the  elections 
still  have  a  finite  energy.  This  is  a  consequence  of 
disregarding  elastic  losses.  Hence  our  analysis  is  not 
applicable  to  very  low  fields  either.  Fortunately,  in  laser 
discharges  the  E/N  values  used  are  such  that  electron  energies 

lie  well  between  these  two  extremes  and  the  approximations 
which  led  to  these  limits  are  acceptable. 

Normalization  provides  one  relation  (18)  between  the 
unknown  excitation  rates  v  and  v_ .  The  detailed  molecular 

X  a 

excitation  processes  considered  in  the  next  section  will  provide 
the  second  relation  whicl  is  necessary  to  determine  them. 


3.  Excitation_of_Diatcmic_Molecul.es 

We  assume  that  the  lasing  molecular  species  can  be 
described  by  a  diatomic,  harmonic  oscillator  model.  The 
Landau  and  Teller7'8  model  gives  the  rate  equation  as  a 
function  of  the  coefficients  for  vibrational-translational 
(V-T)  and  for  vibrational-vibrational  (V-V)  energy  transfer 

g 

in  molecular  collisions .  These  are  introduced  below  : 
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dvc 

n 

at  ' 


[  (n  +  l)e 


'P'l'oT'(n  +  Dxn+. 


-ti 


°  +  n)x  +  ne”  *  c*  , } 
n  n-i 


+  31* 


v-v 


10  {  (  n+  1)  (1  +a)xn+1  -  [  (n  +  l)a  +  n'i  +  c))x 


n 


(20) 


+  n  a  x^}  +  nePaf(n  +  l)xn+1  f(ua  -Va)  .  ( (n  +  x) 


f(ua)  +  n  £(ua  -  Va))Xn  +  n  £(n,)Vl)  -  R(«n<t+1  -  «n>l> 


Here,  xn  is  the  fractional  occupation  of  tie  nth  vibrational 
level  of  the  diatomic  molecule-  "r,v"T  qives  the  rate  of  V-T 
relaxation;  0O  is  related  to  the  gas  temperature  T0  by 

frw  V 

O  *  kr~  =  kT~ 
o  o 

where  hwa  is  the  energy  separation  of  the  vibrational  levels. 
V-V 

ZP10  9Ives  the  rt'te  of  V-V  relaxation;  a  is  a  measure 
of  the  energy  in  vibrational  levels. 


00 

a  ™  E  nx 
o  n 


oo 


The  tcro  proportional  to  n  P  oivcs  the  direct  excitation 

C  cl 

and  de-excitation  of  vibrations  by  electron  impact.  The  ex¬ 
citation  of  the  vibrational  levels  is  known  to  proceed  via  an 
intermediate  ionic  state.  The  ci ess-section  for  vibrational 
excitation  by  electron  impact  shows  a  series  of  peaks  above 
a  certain  threshold  energy,  with  increasing  energy.  The  peaks 
at  higher  energies  correspond  to  excitation  of  the  higher 
vibrational  levels  of  the  molecule  from  its  ground  state. 


Here  we  do  not  consider  precisely  such  a  process.  We 
allow  only  for  vibrational  excitation  by  the  electrons  in 
single  quantum  jumps.  This  assumption  lumps  the  entire  ex¬ 
citation  cross-section  at  the  sharply  defined  energy  u&.  This 
assumption  is  less  severe  than  may  appear  at  first.  The  only 
result  of  the  analysis  that  finally  matters  is  the  net  exci¬ 
tation  rate  of  the  different  vibrational  levels.  The  way  these 
rates  are  produced  may  be  chosen  rather  arbitrarily,  mathe¬ 
matical  convenience  being  one  of  the  criteria  for  the  choice. 
The  mathematically  most  convenient  model  is  to  assume  that 
the  cross-section  increases  with  increasing  n  as  done  in 
Ecr.  (20)  . 

An  equation  for  the  energy  parameter  a  may  be  obtained 
from  (20)  by  a  simple  manipulation.  We  multiply  each  equa¬ 
tion  in  the  set  by  n  and  sum.  When  the  summations  are 
evaluated,  the  V-V  transfers  disappear  because  they  do  not 
change  the  total  energy  in  vibrational  whatever  their  dis¬ 
tribution.  The  result  is 


a-a 

da  .  ©  v  M  R 

3t  +  "TT  Va  +  R 


(22) 


where  we  have  defined  the  relaxation  rates 


i-  5  SP^U  -  e  °)i  va"  -  nePaf  (1  +  a)  £(«,)  - 


-0 


M 


«£ (u,  -  V,  >  , 


(23) 


and  a  is  defined  by 
o 
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Equation  (22)  is  essentially  the  energy  conservation 
equation  of  molecular  excitation. 


If  we  had  assumed  an  excitation  cross-section  for  vibra¬ 
tional  excitation  by  the  electron  that  is  independent  of  vi¬ 
brational  level,  we  would  i.ave  had  to  make  the  additional 
assumption  of  a  Boltzmann  distribution  of  the  vibrational 
states  in  order  to  carry  out  the  summation.  Upon  eva?.ua'ion , 
we  would  have  found 

v  M  =  — { (1+a)  f  (u  )  -  a  f(ua-Va)> 

a  1+CL 

a 

Hence  ?a  of  the  previous  expression  is  replaced  by  — -  . 


Equation  (22)  then  incorporates  everything  we  need  to  know 
about  the  effect  of  the  electrons  and  the  laser  field  upon 
the  vibrational  system.  We  have  introduced,  in  the  above, 
two  vibrational  excitation  rates,  vM  and  v,M.  The  number  of 
electron  collisions  per  unit  volume  is 


v  n  b  v  ^  n 
a  e  a 

which  must  be  equal  to  the  number  of  exciting  collisions 
v  **  of  the  molecular  species  of  density 


M  ne  ,H  .  ^e  Va 

1  "  N  VVa  *  K  CM 


m) 
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where  CM  ■  NM/N  is  the  mole  fraction.  Now  va  itself  depends 
upon  the  operating  conditions.  This  dependence  is  discussed 
below. 

Prom  (22)  it  is  clear  that  the  rate  of  energy  supply  to 
the  lasing  gas  by  the  electrons  is 

Va  va  “  nevaVa  "  Va  ne  Pa  NM  { (1  +“>f<“a> 
-a£(ua-va)}  (25) 

^ow>  from  Eq .  (16)  we  obtain  f  (u  )  and  f (u  —  V  )  for  aiven 

a  a  a  J 

vaVa.  Introducing  this  into  Eq.  (25)  and  solving  for  v  results 


v 

-a  - 


V2ud 


l+f*aa^  ^x/2ud>  r. 


where 


2u  (1  >  u'  «  *  u  r  =^2. 

a'  Yu  '  a  la  a 


,  3W  ua 

^  a  “  IJ’I  V~  (27) 

2vc  ua  f  x 

The  model  with  excitation  cross-sections  independent 
of  vibrational  level  would  have  resulted  in  a  similar  ex¬ 
pression,  with  £  replaced  by  £  /(1+a). 

a  a 

Containing  the  normalization  (18)  with  (26) ,  one  finds  expressions 


vu2ud  2ufl  +  rayx 
v*  2ud  +  V(va  ^  raVx> 


vc2ud 


2ud  +  WVa  +  or  V  ) 
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Some  interesting  features  may  be  noted  in  expressions 
(28)  and  (29)  .  For  a  very  large  ud  (or  E/N) ,  and  a  =  0, 
v  approaches  £  v_.  Thus  £  expresses  the  rate  of  vibra- 
tional  excitation  at  high  fields  in  terms  of  the  elastic 
collisions  when  all  molecules  are  in  the  ground  vibrational 
state , 

Next,  vc  note  that  v_  and  v  approach  zero  when  E/N  0 

21  X 

(u^  0)  .  This  is  clearly  as  it  must  be.  From  this  fact 

2 

one  may  obtain  a  check  on  (28)  and  (29).  Note  that  vy  *  ud 

and  v  «  u , ,  when  a  =  0  and  u ,  -*■  0  .  Because  of  the  nor- 
a  a'  d 

realization  (18)  ,  we  find 


which  checks  (29)  for  a  >=  0  and  ud  0 . 

It  is  further  of  interest  to  note  that  v  decreases 

21 

with  increasing  a,  for  fixed  driving  field  fixed  u^.  This 
is  easily  understood  from  the  fact  that  superelastic  collisions 
come  into  play  when  a  /  0.  These  collisions  impart  energy  to 
the  electrons  at  the  expense  of  the  vibrational  mode  thus 
decreasing  the  net  rate  of  energy  supply  by  the  electrons, 

w 

'*x,  on  the  other  hand,  is  enhanced  by  increasing  a, 
particularly  at  low  values  of  ud.  hhen  ud  is  low,  f (u) 
has  a  larger  discontinuity  at  u  =  ufl  than  when  ud  is  high. 

An  increase  in  a  decreases  the  discontinuity  and  "lets  more 
electron  through  to  the  energy  uv,  at  which  electronic 
excitations  occur,  hence  increasing  the  rate  v  . 


4.  Comparison  Wx  th_Niqhan 'sRe suits 

Nighan  has  presented  detailed  ccmoutations  for  the 
electron  distribution  in  CO.  We  ^/ant  to  study  how  the 
analysis  presented  here  meshes  with  Nighan’s  results.  In 
order  to  make  approximations  equivalent  to  his,  wa  assume 
that  all  molecules  are  in  the  ground  state,  a  =  0.  From 
(28)  and  (29) ,  we  obtain: 

2Uj 


V 

vx  =  vc  — *. 


1+c  va 
1+5a  Ya 


(30) 


va  ■  *avc  ,  -  va 

1+ia  Yq 


The  excitation  rate  at  large  E/N  (r)r  ..  .  .  £  \>  ...  .  , 

or  Ud  1S  a  c*NlCfhan  s 

Fig.  12  gives 

v  I  '  3  *  10'8  Nh  - 
E-h» 

VJe  may  use  this  one  result  to  find  all  other  dependencies. 
Thus,  for  example,  the  fractional  power  transfers 


v  V 
a  a 

v  v  +v  V 


and 


xx  a  a 


v  V 

XX  , 

1 

xx  a  a 


v  V 
a  a 


v  V  +v  V 
xx  a  a 


(31) 


behave  qualitatively  like  those  shown  in  Fig.  8  of  Nighan1  for 
pure  CO.  The  two  are  equal  at  a  value  of  u^  aivon  bv 


2ud  "  *aVa 
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Introducing  the  definition  o'  ud  one  has  for  the  field  at 
which  the  two  curves  cross 


(4, 

N' cross 


Vc 

N 


V 

a 


(32) 


We  use  vc/N  *  3.2  x  10-7  (the  cross  section  for  elastic 
collisions  in  ~0  peaks  at  Pc  =  120  at  an  electron  energy 
1.7  eV  according  to  reference  9,  and  hence  £a  «  3/32.  The 
collision  freouency  Is  v  _/N  =  v  P_p/N,  where  p  is  the  pres- 

C  V 

-7 

sure  in  torr.  h  peak  collision  rate  of  3.2  x  10  results ,) 
Further  V.  *  .2  volts.  Thus, 

a 


(E) 

N  cross 


-  1.06  x  10"15 


volt  cm 


2 


which  compares  well  with  Nighan's  results.  Figure  5  shows 
the  plot  of  the  fractional  power  transfer  to  molecular  excita¬ 
tion  and  Nighan's  curve  for  comparison.  The  agreement  is 
reasonable;  the  E/N  values  are  off  by  about  20%  and  the  effi¬ 
ciency  does  not  decrease  as  fast  in  the  simplified  model. 

This  latter  effect  may  be  understood  on  the  basis  of  the  fact 
that  the  simplified  model  possesses  a  "sink"  in  energy 
space  at  u  -•*  u  and  hence  forces  the  electron  distribution 
to  remain  at  lower  energies. 

Next  consider  the  plot  of  Veff/Nj  of  Nighan's  Fig.  12 
and  the  plot  of  vfl/CMN  obtained  from  the  simplified  model. 
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Figure  6  shows  such  a  comparison.  Again  we  find  qualitative 
agreement.  The  energy  range  of  the  electrons  is  restricted 
to  2ufl/5  <  u  <  ux/5,  which  accounts  for  the  sharper  dropoff 
of  veff/N  at  low  energy  levels  and  the  termination  of  the 
curve  at  uft/5  as  obtained  from  the  simplified  model. 

5.  Electron_Distribution_and_Pum£in2_in_a_C02  Laser  Model 

In  sections  2  and  3  the  electron  distribution  in  a 
diatomic  gas  laser  was  obtained  and  its  interaction  with 
the  population  distribution  of  the  lasing  species  was 
determined.  In  the  C02  laser,  the  lower  laser  level  per¬ 
tains  to  a  different  vibrational  mode.  In  order  to 
analyze  laser  operation  ".11  modes  of  vibration  have  to 
be  taken  into  account  separately.  In  this  section  we 
extend  the  analysis  in  this  way . 

The  model  we  adopt  here  for  mathematical  simplicity  is 
one  in  which  the  molecular  gas  is  represented  by  two  tempera¬ 
tures.  For  this  purpose  we  consider  the  nitrogen  vibrational 
levels  tightly  coupled  to  the  asymmetric  stretching  mode  of 
C02  and  treat  the  two  vibrational  systems  as  one.  In  a  simi¬ 
lar  way,  we  assume  that  the  symmetric  stretching  mode  of  Z02 
is  tightly  coupled  to  the  bending  mode  and  that  both  of  these 
are  describable  by  a  single  vibrational  temperature.  We  fur¬ 
ther  assume  that  the  nitrogen  vibrational  mode  and  the  asym¬ 
metric  stretching  mode  of  CC>2  are  excited  via  an  intermediate 
ionic  state  so  that  excitation  occurs  when  the  electrons  have 


reached  a  critical  energy  ufl  (=2eV);  in  excitation  they 
are  assumed  to  lose  a  portion  of  this  energy.  The  ex¬ 
citation  of  the  combined  symmetric  stretching  mode  and 
bending  mode  is  assumed  to  occur  by  low  energy  electrons 
of  energy  ub  which  upon  excitation  lose  an  energy  Vb. 

When  these  simplifying  assumptions  are  made,  closed  form 
expressions  can  be  derived  for  all  physical  quantities  of 
interest. 

The  distribution  function  of  the  electrons,  f(u), 
again  obeys  Eq.  (15)  .  We  supplement  the  G  function  re¬ 
presenting  the  effect  of  inelastic  collisions  by  the  con¬ 
tribution  of  the  combined  bending  symmetric  stretching  mode 
excitation  characterized  by  an  excitation  frequency 
The  new  G  function  on  the  right  hand  side  of  (15)  is  shown 
in  Fig.  7 .  In  addition  to  supplementing  G  +  G  by  a 
rectangular  function  Gb  representing  the  bending  mode  colli¬ 
sions,  we  return  the  electrons  which  excite  electronic  levels 
at  u- ,  not  at  u  ■  0 . 

ii  (i;)2vc  “V2  a§  ■  -<V  Ga  +  Gb>  <”>' 

Here  vx  is  the  frequency  of  electronic  excitations  assumed 

+ 

to  occur  when  the  electrons  reach  energy  u:  v_  is  the  exci- 
tation  frequency  of  the  intermediate  ionic  state  of  N2  which 
then  leads  to  an  excitation  of  the  combined  nitrogen  vibra¬ 
tional  levels  and  asymmetric  stretching  mode;  vb  is  the 
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frequency  of  excitation  of  the  combined  bending  and  symmetric 
stretching  modes  which  occurs  when  the  electrons  reach 
energy  u^.  See  Figure  7  for  the  G-f unction.  Introducing 
again  the  definition  for  drift  energy  u^,(17)f  the  closed 
form  solution  for  the  distribution  function  is 


v  3  i  i 

f(u)  =  —  -  (“  -  -=r)  E  f,(u);  u  <  u  <  u* 

Vc  2ud  ^  ^x 


va  3  V 

f  (u)  =  f .  (u)  +  -1 - It*.  =  f ,  (u)  ;  Uv  <  u  <  Ua 

1  v  >2u .  2u  3/2  ~  X  a 

c  a  a 


f(u)  =  f2(ux)=  f3(u); 


ub  <  u  <  ux 


v.  3  1  1  + 

f  (u)  «  f ,  (u)  + - -  -=p)  =  f.  (u) ;  u.  <  U  <  u. 

3  vc2ud  *  ^  4  b  b 


f  (u)  »  f4  (ub) ;  0  <  u  <  ufa 


(34) 


If  one  introduces  the  normalization  condition 


/  f  (u)  /u  du  =  1 
0 

one  obtains  a  relationship  among  the  different  excitation 
frequencies 


Vxvx  +  vava  +  Vbvb  =  vc  2ud 


where 

Vx  -  ux 

and 

vb  “  "b 

(35) 


Equation  (35)  is  a  form  of  the  energy  conservation  law.  On 
the  left  hand  side  appears  the  rate  of  energy  transfer  by 
the  electrons  to  electronic  excitations  and  vibrational  mode 
excitations  and  on  the  right  hand  side  is  the  rate  of  energy 
supply  by  the  electric  field  to  the  electrons.  Equation  (35) 
•  is  ore  equation  for  the  three  unknowns,  v  ,  v  ,  vK  for  a 
given  electric  field  E.  In  order  to  find  more  relationships 
one  has  to  study  the  molecular  excitation  in  greater  detail. 


6.  The  Molecular  Excitation  of  the  C02-N2  Model 

We  assume  that  the  distribution  of  C02  over  the  asymmetric 
stretching  mode  (and  the  nitrogen  vibrational  distribution) , 
the  symmetric  stretching  mode,  and  the  bending  mode  are  all 
described  in  terms  of  vibrational  temperatures.  We  introduce 
convenient  energy  parameters  by 
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b  =  exp  (- 


(36) 


a 


exp(- 


s  =  exp  ( 


) 


where  the  energy  spacing  of  the  harmonic  oscillator 

m^xie  i,  i  *  a,  b,  s. 

The  energy  parameter  for  the  symmetric  stretching  mode, 

2 

s,  is  equal  to  the  square  of  that  of  the  bending  mode,  s  *  b  , 
because  the  assumption  is  made  that  the  spacing  of  the  energy 
levels  of  the  symmetric  stretching  mode  is  twice  that  of  the 
bending  mode  and  that  the  temperatures  of  the  two  modes  are 
the  same.  It  is  helpful  to  introduce  additional  energy  para¬ 
meters  by  the  definition 


a  _ 

i+h:  = 


a 


(37> 


The  energy  U  in  a  particular  mode  is  given  directly  in  terms 
of  the  Greek  letters,  e.g.  =  ava  The  energy  conserva¬ 
tion  relation  for  the  asymmetric  stretching  mode  can  be 
written  particularly  conveniently  in  the  form 


n  v 
e  a 


cx-a  CM 

NCMN(c£  +  +c~  R> 

3  MN 


(38) 


Here  t  is  a  phenomenological  relaxation  time  for  the  asymmetric 
stretching  mode;  R  is  the  rate  at  which  vibrational  quanta  are 
lost  from  the  system  via  radiative  transitions;  ng  is  the  electron 
density;  N  is  the  particle  density;  CMN  the  mole  fraction  of 
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nitrogen  and  COj  because  a  describes  the  combined  energy  of 
the  asymmetric  stretching  mode  and  the  nitrogen  vibrational 
mode.  In  order  to  write  an  energy  conservation  relation 
for  the  symmetric  stretching  mode  and  the  bending  mode  it 
is  necessary  to  take  note  of  the  degeneracy  of  the  bending 
mode.  The  m-th  level  of  the  bending  mode  has  the  degeneracy 
m+1.  Under  these  conditions  the  energy  in  the  bending  mode 
is  2Bhwb.  The  energy  in  the  symmetric  stretching  mode  is 
ohus  ■  2  The  energy  conservation  relation  for  the  com¬ 

bined  symmetric  stretching  mode  and  bending  mode  is 


a  -  a. 


nevb  *  2ncm  < —  (o  +  B)  +  - 

dt  T 


(B  -  BJ 


-  R>  (39) 


Here  we  have  introduced  the  two  phenomenological  relaxation 
times  and  for  the  stretching  and  bending  modes  respec¬ 
tively.  CM  is  the  mole  fraction  of  C02 .  The  arrival  of 
quanta  via  radiative  transitions  is  contained  in  R.  Eauations 
(38)  and  (39),  using  the  fact  that  a  is  related  to  B  by 

T?o  =  (TTb)  2  (40) 

provide  two  relations  for  the  two  unknowns  a  and  B  in  terms 
of  va  and  v^.  Now,  the  excitation  rates  themselves  depend  on 
the  molecular  excitation.  One  of  these  deoendencies  can  be 
taken  directly  from  (26)  .  Indeed,  the  nature  of  the  excitations 
and  de-excitations  at  energies  lower  than  u  does  not  affect 

cl 
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t 


the  solution  for  the  distribution  function  f (u)  in  the 
range  u  >  u  #  and  hence  all  relations  derived  from  this 

cl 

part  of  the  distribution  function  remain  valid.  We  have 
from  Eq.  (26) 


v 


1+5aa 


r 

a 


Now,  consider  the  de-excitation  of  the  (m+1)  th  level 
to  the  m-th  level  of  the  bending  mode  by  an  electron  collision 
as  the  result  of  which  the  m-th  level's  population  is  increased. 

Assuming  that  the  collision  cross-section  is  Q*°(v),  one 
has,  for  the  rate  of  excitation  of  the  fractional  population 
xm  per  de9enerate  mode,  x^g^  (g^  is  the  degeneracy)  by  entry 
of  particles  from  the  level  m+1, 

’m+1  2  in 

(m+1)  gmU  /4tt  *ne  vQ^U(v)  f(v)  dv 

The  factor  (m+1)  arises  from  the  fact  that  we  assume  inter¬ 
action  of  the~.electron  with  the  molecule  via  a  dipole 
moment  of  the  molecule.  Perturbation  analysis  shows  that 
the  cross-section  for  the  m+1  m  transition  is  proportional 
to  m+1.  The  rate  is  further  proportional  to  the  fractional 
population  per  degenerate  energy  level  of  the  (m+l)-th  level. 

The  term  vQ^(v)ne  is  the  collision  freouency  for  a  1  •+  0 
transition.  This  expression  is  conveniently  written  in  terms 
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of  the  distribution  '"nction  f(u) 

ne(r+1)  5^  /  f<u>  Rin<u>  <*u  (42) 

gm+l  10 

where  R1Q(u)  is  /uvqJ°(v)  as  a  function  of  u. 

In  this  format  it  is  easy  to  take  account  of  the  Kle\n- 
Rossland  relation  (detailed  balance)  which  states  that  the  rate 
of  up-transit  ions  from  the  m-th  level  to  the  (m+l)-th  level 
is  related  to  the  rate  of  down  transitions  in  that  the  f(u) 
of  (42>  aas  to  be  replaced  by  f(u  +  Vb)  where  Vb  is  the  energy 
loss  in  an  exciting  collision.  Taking  these  relations  into 
account,  one  has 

• 

^electrons  "  ne(  ^  f  ^  Rio  ^  du 

v  •  0  “ 

m 

*  7“  I(m+1)  /  f  (u  +  Vb)  R10(u)  du  +  m  /  f(u)  R  Q  (u)  du] 

gm  0  0 

0 

«D 

X  , 

+  m  - -  /  f(u  4  Vb)  R10(u)  du) 

m-i  Q 

If  one  assumes  that  the  V-V  coupling  is  so  strong  that  thermal 
equilibrium  is  established  ir.  the  bending  mode,  then  one  may 
set 

xm  -  (m  +  1)  (1  -  b)2  bm 
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and  one  may  obtain  the  rate  of  energy  transfer  from  the 
electrons  to  the  vibrational  mode  by  multiplying  both  sides 
of  (43)  by  gm  and  adding  over  all  m.  The  result  may  be 
identified  with  the  rate  of  energy  supply  per  molecule,  i  .e .  with 
Vbne  vb'^N  CM*  0ne  has'  after  carrying  out  the  summation: 

nevb  * 

iTcJJ  ”  ne  2(I^S)  £  ?  f(u  +  Vfa)  R10(u)  du 

0 

“  b  /  f  (u)  R10(u)  du]  (45) 

0 


As  a  simplifying  assumption  one  assumes  that  R10 (u)  is  a 
unit  impulse  (delta)  function  at  ufa“  (compare  Figure  5)  . 

R10  *  R  5  £u  ”  ufa“) 


Now,  we  may  introduce  the  expressions  for  f(u)  and  obtain 
for  the  integrals 


/  f(u  +  vfe) 


R10 (u)  du 


_  P  _ 3  V  V  -  V  V  V  - 

'5u ,  v  -  [,_2£  (i  -  _'.x.  )  4  J*  a  x  , 

d  x  V  v  +  '  *  v  Ju~  v  J 

^  5*  C  cl  cl 


(46) 


Here  we  have  introduced  for  convenience  the  normalized  velocities 
vb  •  vx+,etc*  Vlhich  3X6  601:21  to  the  souare  roots  of  the  respective 
energies.  The  otuer  integral  giving  the  superelastic  collision 


rate  is: 


/  f  (u)  R10lu) 


du  ■  R 


3 


?Vx- 


va  va  Vx“ 


(47) 


One  interesting  fact  emerges  from  the  a.'.ove  relations.  The 

vx”  vx~ 

superelastic  collision  rate  contains  the  factor  vb+ 

in  the  square  brackets  that  can  be  made  much  larger  than  ui  ity 

4*  - 

if  vb  “vb  vx“  »  1.  This  is  the  case  when  the  electrons 

experiencingV a” vibrational  collision  are  returned  at  energies 

much  smaller  _than  those  experiencing  an  electronic  collision, 

and  if  Vk  - -  is  not  much  less  than  unity,i.e.  the  rectangle 
vb+ 

representing  G.  is  wider  than  the  distance  of  its  left  hand 
edge  from  the  origin  when  plotted  against  Ai. 

Introducing  these  expressions  into  (45)  one  hus  for 


(v  +v  r  )V 
'vx  a  ax  x 


vb  "  *b(1  '  b  1  li-b)2ud+(i+b)b,brb;i^ 


(48) 


where  we  have  defined  the  parameters 
6NC, 


'a 


~  f/Vx  ~Vx 
<Vc  \  v+  v" 

'  X  X 


vx  vx 


ax 


2ua  Vvx  ”  vx) 


bx 


vb :  vb 

vb  vb 


v+  vx 

X  x 

v+  -  v” 

X  X 


(49) 


The  meaning  of  Kh  nay  be  gleaned  from  some  simple  obser¬ 
vations.  Note  that  (29)  shows  that  in  the  limit  of  large  fields 

(u.. 


^avc 


+  ®)  and  zero  vibrational  excitation  (ct  0)  ,\>a 


~2" 


Now,  in  the  limit  of  large  fields,  in  the  absence  of  vibra¬ 
tional  excitation  (a  =  0,  b  «  0)  ,  (48)  shews  that 


vbl  u 


cU-oo 


c^b 


(50) 


The  meaning  of  is  therefore  analogous  to  that  of  ?a;  it 
expresses  the  ratio  v^/v  at  high  fields,  in  the  absence  of 
vibrational  excitation.  Equations  (26),  (35),  and  (4  8)  give 
three  relations  for  the  three  unknowns  vx>>a  and  v^.  For 
convenience  we  repeat  them  below  ir.  modified  form 


V  v  +  V  v  +  V.  v.  <=  2u  ,v 
xx  aa  d  r>  a  c 


-  Vx  ♦  [- 


2u. 


+  ?aaralva 


(51) 


2u,  1  b 

"5bvx  ~  5braxva  +  l~  ~  +  Kb  ~  rbxlvb 


where 


<  <  -  !a_ 

ra  =  v-  vj  -  v;  ' 

A  certain  systematic  variation  may  be  observed  is  equations  (51) . 
The  diagonal  coefficients  contain  a  and  b  as  "feedback",  in¬ 
creasing  the  multipliers,  and  hence  decreasing  the  corresponding 
rates,  as  the  vibrational  populations  increase.  The  effects  on 
the  stretching  mode  and  bending  mode  are  different,  because  dif¬ 
ferent  deqeneracies  have  been  assumed  for  the  two.  The  "feed- 

2u. 


back"  is  smaller  the  larger 


i.e.,  the  larger  the 
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electric  field.  Solving  for  vx#  va,  and  vb  one  finds 


1  + 


b 

Ju 


db 


[ 


1  + 


C  r  V 
ax  x 


Tu 


da 


:1 


2u, 


Zv*  vk  uda 


Vb  Vx 


i  j.  ~*da  .  **b  b  —  .  -  - 

1  +  or +  r  ^  as  +  c»ra*  vr  3u 


db 


1  + 


<2uda  +  SaV 


(-■  2) 


where  2uda  =  2ud  +  *aaraVx 

and  2udb  =  2ud(TiE)  ♦  5b  x=E  rbxVx 


(53) 


As  2udvc  represents  the  power  absorbed  per  electron  these 
ratios  are  reciprocals  of  the  partition  of  power  between 
electronic,  asymmetric  mode  and  bending  mode  excitations. 
No  other  power  loss  has  been  included  in  the  theory. 


It  is  worth  studying  some  properties  of  the  rates  in 
the  limits  of  large  and  small  E— field.  Consider  first  the 
limit  of  large  E,  ud  -  «  and  similarly  u^--  udb  « 


lim  % 


L-  ^db 

*b  ud 


When  the  levels  are  unpopulated,  udb  *»  ud  and  we  find  that 


lim  vb  -  f,bvc;  b  *  0 
E-+«° 

which  checks  with  Eq-<50)  .  then  the  occupation  of  the 

levels  is  nonzero,  udb  <  ud  and  hence 
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the  rate  vb  is  larger  than  the  vfa  for  b  »  0.  This  is  the 
consequence  of  the  fact  that  the  average  cross  section 
increases  with  increasing  vibrational  temperature.  Even  >  ' 

though  the  superelastic  collision  rate  goes  up  also 
large-field  limit  the  inelastic  collision  contribution  dominates 
Consider  next  v_  in  the  same  limit 

a 


c 

lim(\T) 

E-*-«  a 


1  2u 

1  da 


When  the  levels  are  unpopulated, 


u 


da 


and 


lim  v 
E-*«  a 


hvc 


(55) 


analogous  to  {54)  .  When  the  levels  are  populated,  u,  >  u 

da  d 

and  the  rate  \>a  goes  down.  Since  we  have  assumed  that  the 
levels  associated  with  are  nondegenerate,  the  collision 
cross  section  does  not  increase  as  rapidly  as  that  of  the 
bending  mode,  and  superelastic  collisions  are  relatively  more 
important.  They  are  responsible  for  the  decrease  of 
increasing  a . 

Finally,  the  limit  for  v  is  simply 
2Uj 


va  with 


Next  consider  the  limits  of  the  same  rates  as  E  -►  ud  ■*  0 

Here,  one  has  to  decide  on  the  limiting  behavior  of  uda  and 

u  as  u.  +  0.  This  behavior  follows  from  the  recognition, 
db  d 

confirmed  by  the  evaluation  to  follow,  that  va  is  propor¬ 
tional  to  u  2,  whereas  v  is  proportional  to  ud.  Therefore, 

a  2 
a,  which  is  proportional  to  v&,  is  proportional  to  ud  and, 

therefore,  uda  as  ud  ■►0.  Using  this  fact,  one  finds 


lim  vt 


lim  v. 


lim  v„  « 


4uabud 

VaxVx 


8ud  uda  “db 
«a«b  VbVx  rax 
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7 .  The_V-I_Characteristic  _for _an  J>Bear_Laser  '  11#  12  ^ 


In  this  section  we  shall  apply  the  previously  developed 
formalism  to  an  E-beam  laser.  Because  certain  approximations 
have  already  been  made  in  obtaining  Eqs .  (52)  ,  the  present 
analysis  will  apply  only  to  the  quasi-CW  operation  of  the 
laser  in  which  the  two  temperature  description  of  the  vibra¬ 
tional  modes  is  adequate.  Clearly,  an  extension  of  the  ana¬ 
lysis  presented  here  could  cover  the  five  temperature  model 
with  attendant  complications .  Also,  we  have  damped  the 
gas  temperature  to  room  temperature . 

The  electron  density  in  an  E-beam  laser  is  prespecified 
by  the  E-beam  excitation.  This  makes  the  determination  of 
the  V-I  characteristic’ particularly  easy. 

From  power  conservation  one  has 


E  J  -  e  ne  <Vxvx  +  Vava  +  Vb)  (56) 


which  serves  to  evaluate  J  in  terms  of  the  rates  determined 
earlier.  We  introduce  the  normalization  condition  (35) 
and  obtain 


J 


2u.v„en^ 
a  c  e 


(57) 


7 

Now,  nfi  is  prespecified,  u^  «  E  .  Therefore,  the  relation 
between  J  and  E  is  a  linear  relation  with  a  given  by 


eV 


mv 
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regardless  whether  lasing  occurs  or  not.  This  fact  is  rather 
important.  We  conclude  that  the  V-I  characteristics  of  the  sustains 
voltage  electrodes  would  remain  unchanged  in  the  presence  or 
absence  of  lasing.  What  happens  is,  that  for  a  given  power 
input  per  unit  volume  all  the  lasing  does  is  to  redistribute 
the  flow  of  energy  among  the  levels,  changes  vx,  va  and 
changes  the  electron  distribution,  but  does  not  affect  the 
conductivity  of  the  medium. 

The  normalization  condition  would  contain  one  additional 
elastic-collision  loss  term.  The  J-E  relation,  however,  would 
still  be  linear  and  have  the  same  form  as  (54) . 

We  may  ask  what  assumption  Jn  the  above  model  led  to 
the  result  of  a  c  independent  of  laser  action.  In  general 
the  conductivity  is  given  by 


a 


2 

3 


nee 


m 


/ 

f  (o) 


2 


df  =  - 
3 


u 


3/2 


f  d 


v 

c 


so  that  the  above  formula  results  for  a  constant  vc  whatever 
the  distribution  f(u).  Lasinq,  or  any  other  energy  loss, 
usually  decreases  the  mean  energy  of  the  distribution.  As 
v  (u)  is  generally  an  increasinq  function  of  u,  lasino  may 
be  expected  to  increase  the  conductivity  and  the  plasma  may 
take  more  power  from  a  constant  voltage  source. 
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®  •  £ain_and_Ef  f  iciencv_qf_an_E-Eeani_Laser 

The  operation  of  an  E-beam  laser  [10 ,  12]  is  modelled 
by  Eqs.  (52)  which  give  the  excitation  rates  v_  and  v 
in  terms  of  the  sustainer  field  E  and  the  temperatures  of 
the  vibrational  modes.  The  conservation  equations  for  the 
molecular  systems (38)  anc  (39)  give  the  vibrational  temperature  fac¬ 
tors  a,b  in  terms  of  vb,  vfi  and  R  —  the  rate  of  lasing 
transitions.  The  gain  constant  a  of  the  P(J) -transition 
of  the  laser  is  obtained  from 

X2  2B 

ag  *  “  AT2  <2j  “  1)  {Nqo1  exp[-BJ  (J  -  l)/kTQ] 

o 

-  N100  exp[-BJ (J  +  1)/VT0]> 

where  A  is  weakly  J-dependent.  Stacz  et  al.  [12]  have 
computed  A  -  .23  sec”1.  Until  now  we  did  not  need  to 
make  any  assumption  about  the  distribution  of  the  mole¬ 
cules  over  the  \ ibrational  states .  The  equations  in¬ 
volved  energy  supplied  or  taken  out  of  the  vibrational 
modes  and  gave  full  information  concerning  these  energy 
rates.  In  order  to  evaluate  the  laser  population  inver¬ 
sion,  it  is  necessary  to  make  the  assumption  that  the 
molecules  are  distributed  over  their  vibrational  states 
according  to  the  Boltzmann  distributions.  Then,  nqoi' 
the  population  of  the  quantum  state  001,  has  the  probability 

(1  -  a)  (1  -  b)2(l  -  s) a 
The  probability  of  the  100  state  is 

(1  -  a)  (1  -  b)2(l  -  s ) s 

2 

s  «■  b  ,  and  thus , 
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N001  ■  N  CM  (1  -  a) (1  -  fc)2 (1  -  b2) a 


and 


N100  *  N  cm  (1  "  a)(l  -  b)2(l  -  b2)b2 

We  may  thus  compute  the  gain  as  a  function  of  the  electric  field 
and  R.  We  show  the  result  in  Fig.  8  for  no  laser  radiation, 

R  »  0.  The  values  of  the  physical  parameters  used  to  ob¬ 
tain  the  gain  are  displayed  in  Table  I.  The  efficiency  is 

„  N  CM  R _ 

°e(vbVb  +  vaVa  +  vxVxJ 

The  plot  of  n  vs  E/N  is  also  shown  in  Figure  9  for  different 
values  of  percent  ionization  (n^/  N  CM)  for  the  case  when 
R  is  such  as  to  halve  the  gain,  and  also  when  R  is  so 
large  as  to  have  reduced  the  gain  to  zero. 

Figure  10  shows  the  rates  vh,  va  and  vy  as  functions  of 
E/N,  for  varying  percent  ionizations,  for  the  three  values 
of  R:  R  «  0,  R  such  that  cxQ  is  reduced  to  half  its  R  =  0 
value,  R  at  aQ  =  0.  We  note  from  Figures  10  that  the  v 's 
arc  roughly  independent  of  the  percent  ionization  and  the 
rates  of  excitation  per  unit  volume  are  proportional  to  n 
The  rates  "droop"  at  the  lower  values  of  E/N  with  increasing 
laser  excitation,  the  droop  being  the  less  the  larger  the 
Percent  ionization.  This  may  be  understood  from  the  fact  that 
the  "loading"  of  the  electron  distribution  by  the  molecular 

pumping  and  lasing  is  the  less  severe  the  higher  the  percent 
ionization. 
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9  .  Ex£ension_to_Study_of_Se if -s us tained__Dis charge 
Char acteris tics 


The  present  formalism  can  be  extended  to  study  laser 
action  in  self-sustained  discharges.  The  rate  of  ioniza¬ 
tion  is  related  to  the  rate  of  electronic  excitation 
-Bp/E 


vi 


■  kv  e 
x 


(58) 


where  B  is  now  one  of  the  Stoletow  constants,  and  p  is  the 
pressure.  The  ionization  frecp&ncy  itself  is  determined 
by  diffusion  and  recombination  losses: 


Vjl  -  a  n  +  D/A2  (59) 

where  D  is  the  diffusion  constant  or  is  the  recombination 
coefficient  and  A  is  a  parameter  proportional  to  the  dis¬ 
charge  tube  radius  [12]  .  The  other  pertinent  equations  are 
Eqs.  (52)  for  \>fa,  and  vc,  and  Eqs.  (38)  and  (39)  in¬ 
volving  a  and  8 . 


For  a  given  E/N,  laser  radiation  (R)  ,  gas  temperature 

(o  ,  8  ),  and  tube  geometry,  Eq3  .  (52),  (38)  and  (39)  give 

the  five  unknowns  a,  8,  v^,  vft,  and  vx  in  terms  of  nQ.  Eg. 

(58)  gives  then  ne  can  be  detsrmined  from  (59).  If  the 

2 

system  is  diffusion  dominated  [a _n  <<  D/A  )  then  Eq.  (59) 

'  *•  G 


serves  to  determine  the  tube  radius  R  in  terms  of  the  elec¬ 
tric  field  (i.e.  in  this  limit  the  E  field  is  a  unique 
function  of  R  and  D)  .  If  we  want  to  determine  the  current, 

ve  need  an  additional  equation  for  J  which  is  provided  by 
(57). 


10 .  Conclusions 

The  purpose  of  this  paper  was  to  obtain  equations  for 
the  pumping  of  molecular  systems  by  electrons  in  a  discharge, 
and  to  do  this  self-consistent.ly  so  as  to  take  the  reaction 

of  the  molecular  distribution  upon  the  electron  distribution 
into  account. 

It  was  found  that  the  normalization  condition,  which 
is  the  power  conservation  relation  for  the  electrons  driven 
by  the  electric  field,  gave  one  relation  between  the  mole¬ 
cular  excitation  rate(s)  and  the  electronic  excitation  rate. 
The  necessary  additional  relations  for  the  molecular  exci¬ 
tation  rates  were  obtained  by  a  detailed  study  of  the  in¬ 
elastic  and  superelastic  collision  processes.  In  this  way 
relations  were  obtained  for  the  rates  in  terms  of  the  vibra¬ 
tional  population  distributions  (molecular  temperatures)  . 
T.iese  in  turn  must  satisfy  the  energy  conservation  relations 
for  the  molecular  pumping  and  relaxation  processes.  In  this 
way  as  many  equations  as  unknowns  are  obtained  for  the 
analysis  of  a  discharge  with  prespecified  ionization,  such 
as  in  the  E-beam  laser. 
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The  E-beam  laser  provides  a  simple  application  for 
the  formalism.  The  results  obtained  from  the  present  analysis 
are  eminently  reasonable.  They  predict  the  behavior  of  a 
COj  E-beam  laser  qualitatively  and  sven  semi -quantitatively. 

It  la  hoped  that  the  formalism  developed  here  is  simple 
enough  so  that  it  will  lend  itself  to  the  analysis  of  prob¬ 
lems  involving  spatial  boundary  conditions,  flow,  diffusion, 
situations  for  which  a  simplified  pumping  analysis  is  neces¬ 
sary  in  order  to  reduce  the  analytical  task  to  manageable 
proportions . 
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Table  I 


tr  -  10  eV 


Vx  ■  9  eV 


ua  -  2  eV 


Va  -  0.289  eV 


Vb  "  *25  va 


u£  *  0.077  eV 

vc/N  -  0.47  x  10”7  sec"1  cm3  [13] 


Ca  -  0.076 


lAa  *  380  sec 


[11 

-1 


Cb  -  0.107  [1] 


[14] 


1/Tb  *  16,000  sec"1  [14]* 


VV=  0 


BAT  -  0.0019 
o 


!/ir  T,  -  4.74  MHz/torr  [15] 


-  0.12  eV 


He:C02:N2  -  6.6:1: .43 


*  The  two  relaxation  rates  of  the  symmetric  stretching 

mode  and  the  bending  mode,  have  been  replaced  by  a 

single  rate  lAfa,  the  value  of  which  is  ecruated  to  the 

rate  of  relaxation  of  the  symmetric  stretching  mode. 

2 

Since  o*6,  this  corresponds  to  a  linearization  of 
Eo.  (39)  . 
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Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


Figure  Captions 

1  Role  of  s^nei  as  source  of  electrons  in  energy 
space 

2  The  negative  divergence  of  G,  +  G„  balances  the 
(positive)  divergence  of  G£ 

3  The  currents  Gv  and  G 

X  cl 

4  The  distribution  function  for  punping  of 
an  idealized  diatomic  gas 

5  Nighan's  plot  of  efficiency  of  molecule-  excita¬ 
tion  compared  with  present  results 

6  Comparison  of  Nighan's  vef^  with  of  the 
present  analysis 

7  The  current  in  energy  space,  GE  for  two  mole¬ 
cular  excitation  cross-sections 

8  The  small  signal  gain  of  E-beam  laser  with 
percent  ionization  as  parameter  (the  gas  tempe/a- 
ture  is  assumed  to  be  room  temperature) 

9  Efficiency  n  vs  E/N  with  percent  ionization 
as  parameter 

10  The  rates  of  excitation,  vfa,  va,  and  vx 
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Figure  9 
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